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摘 要 :针对 频率 分 集 阵列 三 维 合 成 孔径 雷达 (3D -FDA-SAR) 成 像 方法 阵 元 数目 多 、 阵 元 利用 率 低 的 问题 ,提出 一 种 基于 
多 输入 多 输出 阵列 的 3D-FDA-SAR 成 像 方法 。 将 3D-FDA-SAR 切 航 迹 向 的 频率 分 集 阵 列 改 为 多 输入 多 输出 频率 分 集 阵 
列 , 多 输入 多 输出 频率 分 集 阵列 随 运动 平台 运动 ,在 沿 航 迹 向 形成 综合 孔径 ,结合 切 航 迹 向 的 实 阵列 ,合成 虚拟 二 维 频率 分 
集 阵 列 平面 ,获得 对 目标 的 下 视 三 维 成 像 能 力 。 首 先 建立 多 输入 多 输出 3D -FDA-SAR 成 像 模型 与 信号 模型 ,利用 多 输入 
多 输出 技术 ,在 发 射 端 发 送 波形 正 交 信号 单 频 窄 带 信和 号 ,在 接收 端 通过 全 频 接 收 方式 ,接收 所 有 发 射 阵 元 经 目标 反射 后 的 
回 波 信号 ,经 正 交 匹配 滤波 器 进行 信号 分 离 , 获 得 不 同 收发 通道 的 回 波 数 据 ,再 利用 后 向 投影 算法 对 回 波 数 据 进行 成 像 处 
,最 终 获得 目标 的 三 维 成 像 结果 。 仿 真 结 果 表 明 ,基于 多 输入 多 输出 阵列 的 3D -FDA-SAR 成 像 方法 使 用 少量 阵 元 ,提高 
车 元 利用 率 ,获得 了 对 下 视 目 标的 三 维 成 像 能 
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bstract: Aiming at the problem of the large number of array elements and low utilization rate of the three-dimensional syn- 


ssthetic aperture radar (3D-FDA-SAR) imaging method of frequency diversity arrays. a 3D-FDA-SAR imaging method based 


multiple input and multiple output arrays is proposed. The frequency diversity array of 3D -FDA-SAR tangent track di- 
NS on is changed to multiple input multiple output frequency diversity array. The multiple input multiple output frequency 
diversity array moves with the moving platform to form a synthetic aperture along the track direction, combined with the 
tangent track direction The real array is combined with a virtual two-dimensional frequency diversity array plane to obtain the 
downward-looking three -dimensional imaging capability of the target. Firstly. a multi-input multi-output 3D-FDA-SAR im- 
aging model and signal model are established. using multi-input multi-output technology. the waveform quadrature signal 
single-frequency narrowband signal is sent at the transmitting end, and all the transmitting arrays are received through the 
full-frequency receiving mode at the receiving end. The echo signal reflected by the target is separated by a quadrature 
matched filter to obtain the echo data of different receiving and sending channels. and then the echo data is imaged by the 


backward projection algorithm. and finally the three -dimensional imaging result of the target is obtained. Experimental sim- 


ulation results show that the 3D -FDA-SAR imaging method based on multiple -input multiple -output arrays uses a small 
number of array elements, improves the utilization of the array elements. and obtains the three -dimensional imaging capa- 


bility of downward-looking targets. 
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合成 孔径 雷达 (synthetic aperture radar, 简称 
SAR) 包 括 二 维 SAR 和 三 维 SAR, 三 维 SAR 克服 了 
ZHE SAR 存在 的 车 掩 等 问题 ,成 为 当前 SAR 研究 的 
热点 。 三 维 SAR 成 像 技 术 主 要 包括 干涉 SARH LJ 
lr SAR” 、 圆 周 SAR 、 阵 列 三 维 SAR 等 。 与 其 
他 三 维 SAR 相 比 ,阵列 三 维 SAR 通过 在 运动 平台 切 
航 迹 向 放置 线 阵 , 线 阵 随 平台 在 空间 中 运动 ,形成 沿 
航 迹 向 的 综合 孔径 ,结合 切 航 迹 向 的 实 孔 径 , 合 成 二 
维 虚 拟 阵 列 平 面 ， 再 利用 距离 向 的 脉冲 压缩 技术 获得 
对 目标 的 三 维 高 分 辨 成 像 能 力 , 具 有 简单 灵活 的 特 
点 ,成 为 三 维 SAR 重要 研究 方向 。 
= 传统 阵列 三 维 SAR 基于 相 控 阵 雷 达 体制 ,射频 
前 器 发 送 宽带 人 言 号 ,由 此 带 来 了 宽带 信号 非 线 性 和 接 
改 滑 宽带 信和 号 不 人 EARO ERI RISE, 2006 年 , An- 
tonik AU pa h T dui RAE REI] (frequency diverse 
a 和 简称 FEDA) 概 念 ,发 射 阵 元 不 再 发 送 大 时 宽 辜 
究 辜 号 ,而 是 发 送 具有 频率 差 的 音频 信号 ,其 带宽 为 
发 射 阵 元 不 同 单 频 分 量 拼接 获得 。 同 时 其 波束 方向 
RRAN, EB AERA ARAS ,因此 其 对 目 
hg tr foni. 存在 距离 和 角度 的 模糊 问题 。 针 
对 前 一 问题 , 顾 地 良 等 中 提出 了 一 种 多 次 改变 频率 
分 集 阵 列 频 偏 的 方法 ,并 基于 MUSIC 算法 对 多 目标 
进 每 定位 。 文 献 [12] 采 用 随机 频率 增 量 的 方法 ,获得 
理想 图 钉 状 波束 方向 图 ， 并 使 用 压缩 感知 算法 进行 成 
像 仿真 验证 。 文 献 [13] 进 一 步 对 随机 频率 分 集 增 量 
系数 进行 分 析 ,提出 了 一 种 多 目标 无 模糊 定位 方法 。 
频率 分 集 阵 列 三 维 成 像 是 二 维 成 像 的 自然 延伸 ,文献 
[14-16 利用 二 维 线性 频率 分 集 平 面 , 获 得 对 目标 的 
三 维 成 像 能 力 , 但 其 波束 方 向 图 是 三 维 耦 合 的 。 进 一 
步 的 ,Li 等 玉 提 出 频率 分 集 阵列 三 维 合成 孔径 雷达 
成 像 (three dimensional frequency diverse array syn- 
thetic aperture radar. 简称 3D -FDA-SARO Jr i. 3C 
现 了 对 目标 无 模糊 成 像 , 然 而 该 方法 还 存在 发 射 阵 元 
数目 多 、 阵 元 利用 率 低 的 问题 。 多 输入 多 输出 合成 孔 
TE TRACTU ,发 射 端 发 送 正 交 信和 号 ,接收 端 进行 正 交 
信号 分 离 ,使 用 少量 阵 元 ,可 等 效 多 个 相互 独立 的 虚 
拟 阵 元 ,提高 了 阵 元 利用 率 。 鉴 于 此 ,将 多 输入 多 输 
出 技术 引入 3D-FDA-SAR 系统 中 ,提出 多 输入 多 输 
出 3D-FDA-SAR 成 像 方法 。 


1 多 输入 多 输出 阵列 构 型 
在 系统 几何 构 型 满足 远 场 条 件 下 ,一 对 收发 分 置 


的 天 线 阵 元 ,可 以 由 其 基线 中 心 位 置 的 一 个 虚拟 收发 
共用 的 阵 元 来 等 效 。 多 输入 多 输出 阵列 采用 发 射 阵 
元 放置 在 阵列 两 端 , 接 收 阵 元 放置 在 阵列 中 间 的 阵列 
构 型 。 图 1 以 4 发 4 收 为 例 ,展示 了 在 该 阵列 构 型 
下 ,实际 收发 阵 元 位 置 和 等 效 阵 列 示 意图 ,其 中 发 射 
阵 元 间距 为 24 , 接收 阵 元 间距 为 4d, 收发 阵 元 间距 
与 等 效 阵列 间距 均 为 d。 


实 阵 列 
等 效 阵 列 
e. " e 
发 射 阵 元 接收 阵 元 ”等 效 收 发 共用 阵 元 


1 多 输入 多 输出 阵列 构 型 及 等 效 阵列 分 布 示意 图 


为 获得 均匀 分 布 的 等 效 阵 列 ,该 阵列 构 型 需 
满足 : 
ls — 2d, 
4dg — Pd. (D 
drr =d, 
其 中 : dr 为 发 射 阵 元 间距 ; dr 为 接收 阵 元 间距 ; dar 
为 收发 阵 元 间距 ; P 为 发 射 阵 元 总 数 ; d 为 虚拟 阵列 
的 阵 元 间距 。 实 际 阵列 长 度 


|, —CP c- PQ —2)d. (2) 
其 中 Q 为 接收 阵 元 总 数 。 等 效 阵 列 的 长 度 
lL, —CPQ — Dd. (3) 
阵 元 利用 率 
21-7. (4) 


2 基于 多 输入 多 输出 阵列 3D-FDA-SAR 成 
像 模 型 


图 2 为 基于 上 述 阵 列 的 3D -FDA-SAR 成 像 几 
何 模型 和 重要 参数 。 其 中 多 输入 多 输出 频率 分 集 阵 
列 装载 于 运动 平台 并 沿 y 轴 ( 切 航 迹 向 ) 均 匀 分 布 。 
BIERE- ERE H, W r 轴 ( 沿 航 迹 向 ) 直线 飞 
行 。 在 不 同时 刻 平台 运动 到 沿 航 迹 向 的 不 同 观测 点 ， 
在 每 个 观测 点 为 发 射 阵 元 重新 选择 一 组 无 重复 的 PP 
个 随机 频率 ,其 中 发 射 阵 元 的 随机 频率 增 量 从 增 量 集 
f {0,1,…,(L 一 1)Af) 中 随机 选取 , 工 表示 增 量 系 
数 的 最 大 值 ,满足 (L DAF « foo 


观测 点 2 观测 点 mm 


A 观测 点 1 


运动 方向 
—— 


e. " 
发 射 阵 元 ”接收 阵 元 


图 2 基于 多 输入 多 输出 3D -FDA-SAR 成 像 模型 
多 输入 多 输出 频率 分 集 阵 列 随 运动 平台 在 0 时 


交办 高 度 H 处 以 速度 ， 沿 * 轴 运 动 。 在 沿 航 迹 向 
ERIS Ju di 的 M 个 观测 点 。 在 第 m 个 观测 点 


b 个 发 射 阵 元 的 频率 为 
Sno Sfo DAP ag =fo iua (5) 


YER. p BRL, WRN 


f^ 


HR ES EL, 


和 由 ,1 … 江 一 1}5A7 为 单位 频率 增 量 。 在 第 m 个 
观测 点 ， 第 p 个 发 射 阵 元 发 送 的 信号 为 
eS, = 各 sos P pee qua (0 


Hog: T, 为 脉冲 持续 时 间 ; $,, 为 发 射 归 一 化 包 络 


的 焉 交 波 形 。 不 同 阵 元 发 射 的 信号 彼此 正 交 ,满足 如 
ER 
e ian UE i G r)dt —0,. pi Pa* 


Tp 


CD) 


las GO? G—o0dt—l.pi—ps. 
" 


T 


el 
f 
© 


其 中 : r 为 时 延 ; pipo 为 发 射 阵 元 ; * 表示 共 斩 运 
算 。 考 虑 远 场 情 况 , 则 在 第 个 观测 点 时 ,第 q 个 阵 
元 接收 到 的 回 波 信号 为 


yma (t) 一 

E S hui a E 

k=1 d 0 

es (Doe (600 hut — cy) Ima tpa P) + 
0 

Ema t) t (8) 


其 中 : ok) 为 第 & 个 目标 的 反射 系数 ; re —2r/c 为 
公共 时 延 ; em GO 为 加 性 噪声 rw O0 为 目标 到 发 
射 和 接收 阵 元 的 往返 时 延 ,包括 目标 到 达 发 射 阵 元 造 
成 的 时 延 r,，, (RO 和 目标 到 达 接 收 阵 元 产生 的 时 延 
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其 中 : r? O 为 发 射 阵 元 与 目标 点 的 距离 ; r? (6) 
为 接收 阵 元 与 目标 点 的 距离 ， 
r (k)= 


(P. 0) —md, Y. FP, Rk) 


=y FP. H’ ， 


a0) 

rt - 
(P, —md, Y FP, y) P.O- HY, 
(11) 


其 中 : P,(k)、.P,(k)、P.(k) 为 第 & 个 目标 在 直角 坐 
标 系 下 的 坐标 值 ; yw 、y 为 发 射 与 接收 阵 元 的 y 轴 
坐标 值 。 在 第 m 个 观测 点 ,经 匹配 滤波 $,,0 一 
ro)e ,得 到 此 观测 点 第 q 个 接收 阵 元 ,接收 来 自 
第 p 个 发 射 阵 元 的 回 波 信 号 为 

Yma) = Y (o (b) e mi at? 


k=1 


P^ Bau atis 
(12) 


3 成 像 方法 与 仿真 结果 


3.1 基于 后 向 投影 算法 的 成 像 方法 
根据 式 (12) 得 到 第 m 个 观测 点 处 的 回 波 数 据 : 


Y m,1.1 Y m,1.2 Ym.1.Q 
Vm,2,1 Y m,2,2 Ut Vm,2,Q 

yn 一 . (13) 
m.P.1 Ym,P,2 Ym.P.Q 


对 成 像 区 域 划分 网 格 体 ,第 g 个 网 格 体 与 第 
个 观测 点 的 第 zp 个 发 射 阵 元 、 第 q 个 接收 阵 元 的 距离 
分 别 为 rmp Cg) era (Os 收发 时 延 为 Tapa) 
其 中 ， 


Emi (g) +r æ) 


Tmp E) 一 2 
对 回 波 数据 进行 相位 补偿 ,得 到 第 g 个 网 格 体 


与 第 m 个 观测 点 , p 个 发 射 阵 元 9 个 接收 阵 元 的 成 
像 处 理 值 : 
人 
遍历 所 有 收发 阵 元 ,得 到 第 g 个 网 格 体 与 第 mm 
个 观测 点 所 有 收发 阵 元 的 相位 补偿 矩阵 为 


10 桂 


DD,,(g)= 
en d Afr, 1 de d Afr, 2 er d Afr, ÍQ 
JYnl,, 2 Mt, p es 2t ag e P n2 mg 


E di PAS pa Em .PAfTn, Pp,2 em .PAfrm, P,Q 


(15) 
联 立 式 (13) 与 式 (15), 得 到 第 m 个 观测 点 处 相位 补 
偿 后 的 成 像 处 理 值 : 
IT,(g)=y,(g)° ©D,(g), (16) 
Hop. 0 表示 哈达 玛 积 ,每 个 网 格 体 需 要 进行 M 次 
观测 。 对 每 个 观测 点 相位 补偿 后 的 成 像 处 理 值 进行 
相 参 累加 ,得 到 最 终 成 像 像素 值 : 


NM 一 1 
te=] 5|. (17) 


m-—0 


- 
ni 仿真 结果 


时 息 据 上 述 后 向 投影 算法 思想 及 处 理 流程 ,对 多 输 


sage 1 所 示 。 
e 参数 值 
Conon; 37.5 
DB m TA 
PC 和 速度 /ms ) 10 
ex 轴 方 向 虚拟 阵列 长 度 /m 4 
= 轴 方 向 虚拟 阵列 长 度 /m 4 
(发射 阵 元 数 40 
接收 阵 元 数 36 
最 大 频 偏 系数 750 
单位 频率 增 量 /MHz 1 


场景 设置 为 X 轴 方 向 7.5~9.5 m.Y 轴 方 向 一 1 一 
1m.Z l7 8]—0.5—1.5 m 的 区 域 。 图 2(a) 为 仿真 
的 目标 场景 ,8 个 点 目标 分 别 位 于 1 m 立方 体 顶 点 位 
置 ,假设 目标 散射 系数 为 1,3D -FDA-SAR 与 基于 多 
输入 多 输出 阵列 构 型 的 3D -FDA-SAR 仿真 结果 如 
图 3(b) 与 (c) 所 示 , 其 中 成 像 结果 是 在 成 像 表 面 最 大 
值 的 一 3 dB 处 进行 绘制 。 

图 3 仿真 结果 表明 ,保持 其 他 参数 一 致 , 且 满 足 
3D-FDA-SAR 实 阵列 长 度 与 多 输入 多 输出 3D - 
FDA-SAR 虚拟 阵列 长 度 一 致 条 件 。 多 输入 多 输出 
3D -FDA-SAR 成 像 方法 仅 使 用 76 个 阵 元 , 阵 元 利用 
率 为 98.6%, 可 以 实现 目标 的 三 维 成 像 。 
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(b) 3D-FDA-SAR 


mE m/m 


(c) 基于 多 输入 多 输出 3D-FDA-SAR 


图 3 仿真 场景 与 仿真 结果 


为 进一步 分 析 三 维 成 像 效 果 ,保持 参数 一 致 ,图 
4 为 3D-FDA-SAR 在 目标 点 (8. 一 0.5,0) 处 的 沿 航 
迹 向 - 切 航 迹 向 、 切 航 迹 向 -高 程 向 、 沿 航 迹 向 -高 程 向 
的 二 维 切面 图 。 

从 图 4 可 看 出 ,3D -FDA-SAR 成 像 二 维 切 面 能 
量 主要 聚焦 在 目标 点 附近 ,同时 在 目标 点 附近 旁 办 较 
为 杂乱 。 多 输入 多 输出 3D -FDA-SAR 在 目标 点 (8， 
一 0.5,0) 处 成 像 二 维 切 面 图 如 图 5 所 示 。 

从 图 5 可 看 出 ,多 输入 多 输出 3D -FDA-SAR 成 
像 方法 ,二 维 切 面 图 清晰 且 较 为 规整 ,使 用 少量 阵 元 ， 
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图 4 3D-FDA-SAR 切面 图 


利用 多 输入 多 输出 技术 进行 正 交 信号 分 离 ,提高 了 阵 
元 利用 率 , 可 以 实现 对 目标 的 三 维 成 像 ,与 理论 分 析 
一 致 。 


4 结束 语 


提出 一 种 基于 多 输入 多 输出 阵列 的 3D -FDA- 
SAR 成 像 方法 。 该 方法 将 多 输入 多 技术 引入 3D - 
FDA-SAR 中 ,建立 了 成 像 模型 与 信号 模型 ,提高 了 
阵 元 利用 率 ,实现 了 对 目标 的 三 维 成 像 。 仿 真 结果 表 
明 ,该 方法 在 相同 阵列 长 度 下 ,提高 了 阵 元 利用 率 , 降 
低 了 平台 载荷 ,有 利于 实际 布 阵 与 实际 应 用 。 
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